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Abstract. A recent homogeneous study of outflow activ¬ 
ity in low-mass embedded young stellar objects (YSOs) 
(Bontemps et al. 1996) suggests that mass ejection and 
mass accretion both decline signihcantly with time dur¬ 
ing protostellar evolution. In the present paper, we pro¬ 
pose that this rapid decay of accretion/ejection activity 
is a direct result of the non-singular density profiles char¬ 
acterizing pre-collapse clouds. Submillimeter dust contin¬ 
uum mapping indicates that the radial prohles of pre- 
stellar cores flatten out near their centers, being much 
flatter than p{r) oc r~^ at radii less than a few thousand 
AU (Ward-Thompson et al. 1994). In some cases, sharp 
edges are observed at a finite core radius. Here we show, 
through Lagrangian analytical calculations, that the su¬ 
personic gravitational collapse of pre-stellar cloud cores 
with such centrally peaked, but flattened density profiles 
leads to a transitory phase of energetic accretion imme¬ 
diately following the formation of the central hydrostatic 
protostar. Physically, the collapse occurs in various stages. 
The first stage corresponds to the nearly isothermal, dy¬ 
namical collapse of the pre-stellar flat inner region, which 
ends with the formation of a hnite-mass stellar nucleus. 
This phase is essentially non-existent in the ‘standard’ 
singular model developed by Shu and co-workers. In a 
second stage, the remaining cloud core material accretes 
supersonically onto a non-zero point mass. Because of the 
significant infall velocity field achieved during the first col¬ 
lapse stage, the accretion rate is initially higher than in 
the Shu model. This enhanced accretion persists as long 
as the gravitational pull of the initial point mass remains 
significant. The accretion rate then quickly converges to¬ 
wards the characteristic value ~ /G (where a is the 
sound speed), which is also the constant rate found by 
Shu (1977). If the model pre-stellar core has a Hnite outer 
boundary, there is a terminal decline of the accretion rate 
at late times due to the finite reservoir of mass. 
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We suggest that the initial epoch of vigorous accretion pre¬ 
dicted by our non-singular model coincides with Class 0 
protostars, which would explain their unusually powerful 
jets compared to the more evolved Class I YSOs. We use 
a simple two-component power-law model to fit the di¬ 
agrams of outflow power versus envelope mass observed 
by Bontemps et al. (1996), and suggest that Taurus and 
p Ophiuchi YSOs follow different accretion histories be¬ 
cause of differing initial conditions. While the isolated 
Class I sources of Taurus are relatively well explained by 
the standard Shu model, most of the Class I objects of the 
p Oph cluster may be effectively in their terminal accre¬ 
tion phase. 
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1. Introduction 

Despite recent observational and theoretical progress, the 
initial conditions of star formation and the Hrst phases of 
protostellar collapse remain poorly known. 

It is reasonably well established that low-mass stars 
form from the collapse of centrally-condensed cloud cores 
initially supported against gravity by a combination of 
thermal, magnetic, and turbulent pressures (e.g. Shu et al. 
1987, 1993 for reviews). However, the critical conditions 
beyond which a cloud core becomes unstable and starts 
to collapse are uncertain and still a matter of debate (e.g. 
Shu 1977, Mouschovias 1991, Boss 1995, Whitworth et al. 
1996). In particular, they depend on yet unmeasured fac¬ 
tors such as the strengths of the static and fluctuating 
components of the magnetic field. 

Once fast cloud collapse sets in, the main theoretical fea¬ 
tures of the dynamical evolution that follows have been 
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known since the pioneering work of Larson (1969). Dur¬ 
ing a probably brief first phase, the released gravitational 
energy is freely radiated away and the cloud stays isother¬ 
mal. This initial collapse phase, which tends to produce a 
strong central concentration of matter, ends with the for¬ 
mation of an opaque, hydrostatic stellar object (cf. Larson 
1969; Boss & Yorke 1995). This time is often denoted < = 0 
and referred to as (stellar) ‘core formation’ in the litera¬ 
ture (e.g. Hunter 1977). When the stellar core has fully 
formed, one enters the accretion phase during which the 
central protostar builds up its mass (M*) from a surround¬ 
ing infalling envelope (of mass Mgnv) while progressively 
warming up. The infalling gas is arrested and thermalized 
in an accretion shock at the surface of the stellar core, 
generating an infall luminosity Li„/ « GM*(t)Macc/7?*. 
In the ‘standard’ theory of Shu et al. (1993) which uses 
singular isothermal spheres as initial conditions, the ac¬ 
cretion rate Mace is constant and equal to where 

Oe// is the effective sound speed. The singular isothermal 
sphere, which has p oc is a physically mean¬ 

ingful starting point for the ‘self-initiated’ collapse of an 
isolated cloud core because it represents the (unstable) 
limit of infinite central concentration in equilibrium mod¬ 
els for self-gravitating isothermal spheres (e.g. Shu 1977 - 
hereafter Shu77). Furthermore, Lizano & Shu (1989) have 
shown that magnetically-supported cloud cores undergo¬ 
ing ambipolar diffusion evolve naturally toward a singular 
configuration reminiscent of a singular isothermal sphere 
(see also Ciolek & Mouschovias 1994). However, it is likely 
that actual cloud cores become unstable and start to col¬ 
lapse before reaching the asymptotic singular state, espe¬ 
cially when they are perturbed by an external agent such 
as a shock wave (e.g. Boss 1995, Whitworth et al. 1996). 
In this case, the initial conditions will be effectively non¬ 
singular, and Macc is expected to be time-dependent (e.g. 
Zinnecker & Tscharnuter 1984, Henriksen 1994, Foster & 
Chevalier 1993, McLaughlin & Pudritz 1997). It is this 
possibility that we explore further and compare with rel¬ 
evant observations in the present paper. 

An unsettled related issue concerns the manner in which 
the central hydrostatic stellar core develops. While in some 
models core formation occurs in a dynamical, supersonic 
fashion (e.g. Larson 1969, Foster & Chevalier 1993, present 
paper), this formation is achieved by slow, subsonic evo¬ 
lution of the gas in the scenario advocated by Shu and 
co-workers (see also Ciolek & Mouschovias 1994). An im¬ 
portant observational consequence is that while the dy¬ 
namical models predict the existence of ‘isothermal pro¬ 
tostars’ (in the sense of the initial collapse phase outlined 
above), these do not exist in the standard Shu theory. As 
we will see, in practice both situations probably occur in 
nature. 

Observationally, one distinguishes various empirical 
stages in the evolution of young stellar objects (YSOs) 
from cloud core to (low-mass) main sequence star (e.g. 
Lada 1987, Andre 1994). The youngest observed YSOs are 


the Class 0 sources identified by Andre, Ward-Thompson, 
& Barsony (1993 - hereafter AWB93), which are char¬ 
acterized by very strong emission in the submillimeter 
continuum, virtually no emission below A ~ 10 /im, and 
powerful jet-like outflows. Their very high ratio of sub¬ 
millimeter to bolometric luminosity suggests they have 
Menv » Af*. Thus, Class 0 YSOs are excellent can¬ 
didates for being very young protostars (estimated age 
~ 10"^ yr) in which the hydrostatic core has formed but 
not yet accreted the bulk of its final mass (AWB93). The 
next, still deeply embedded, YSO stage corresponds to the 
Class I sources of Lada (1987), which are detected in the 
near-infrared (A ~ 2 pm) and have only moderate submil¬ 
limeter continuum emission (Andre & Montmerle 1994; 
hereafter AM94). They are interpreted as more evolved 
protostars (typical age ~ 10^ yr) surrounded by both a 
disk and a residual circumstellar envelope of substellar 
mass (~ 0.1-0.3 Mq at most in p Ophiuchi; cf. AM94). Fi¬ 
nally, the most evolved (Class H and Class HI) YSO stages 
correspond to pre-main sequence stars (e.g., T Tauri stars) 
surrounded by a circumstellar disk (optically thick and 
optically thin at A ^ 10 /rm, respectively), but lacking a 
dense circumstellar envelope. 

Many examples of pre-collapse, pre-stellar cloud cores are 
known. In particular, Myers and co-workers have studied 
a large number of ammonia dense cores without IRAS 
sources (e.g. Benson & Myers 1989) which are tradition¬ 
ally associated with (future) sites of isolated low-mass star 
formation (see Myers 1994 for a recent review). These star¬ 
less dense cores, which are gravitationally bound and close 
to virial equilibrium, are believed to be magnetically sup¬ 
ported and to progressively evolve towards higher degrees 
of central concentration through ambipolar diffusion (e.g. 
Mouschovias 1991). 

More compact starless condensations have been identified 
in regions of multiple star formation such as the p Ophiuchi 
main cloud (e.g. Loren, Wootten, & Wilking 1990, Mezger 
et al. 1992b, AWB93). In these regions, the weakness of 
the static magnetic field (Troland et al. 1996) and the 
complexity of its geometry (Goodman & Heiles 1994 and 
references therein) suggest that gravitational forces over¬ 
power magnetic ones and that an external trigger rather 
than ambipolar diffusion is responsible for cloud fragmen¬ 
tation and core formation (e.g. Loren & Wootten 1986). 

In the present paper, we bring together two indepen¬ 
dent sets of observational results recently obtained on the 
density structure of pre-stellar dense cores (Sect. 2.1) and 
on the evolution of protostellar outflows (Sect. 2.2), and 
interpret them by means of a simple analytical theory 
(Sect. 3.3) which sheds light on the accretion histories 
found in the numerical work of Foster & Chevalier (1993 
- hereafter FC93). We compare our theoretical predictions 
with observations in Sect. 4 and we conclude in Sect. 5. 

2. Relevant observations 
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2.1. Constraints on initial conditions: flattened pre-stellar 

density profiles li 689 b 


(b) L15B7 


Recent submillimetre dust continuum mapping shows that 
the radial density profiles of pre-stellar cores are relatively 
steep towards their edges (i.e., sometimes steeper than 
p(r) oc r“^) but flatten out near their centres, becoming 
less steep than p{r) oc (Ward-Thompson et al. 1994 - 
hereafter WSHA; Andre, Ward-Thompson, & Motte 1996 
- hereafter AWM96; Ward-Thompson et al. 1997). A rep¬ 
resentative, well-documented example of an isolated pre- 
stellar core is provided by L1689B, which is located in 
the p Ophiuchi complex but outside the main cloud. In 
this case, the radial density profile approaches p{r) oc 
between ~ 4000 AU and ~ 15000 AU, and is as flat as 
p{r) oc or p{r) oc r~^ ‘^ (depending on the deprojec¬ 

tion hypothesis) at radii less than ~ 4000 AU (AWM96). 
The mass and density of the relatively flat central region 
are estimated to be ~ 0.33 Mq and ~ 2 x 10® cm“®, re¬ 
spectively. 

By contrast, protostellar envelopes are always found to 
be strongly centrally peaked and do not exhibit the inner 
flattening seen in pre-stellar cores. In particular, isolated 
protostellar envelopes have estimated radial density pro¬ 
files which range from plj) oc ® to p(r) oc r~^ over 
more than 0.1 pc in radius (e.g., Ladd et al. 1991; Motte 
et al. 1996). This difference of structure between the pre¬ 
collapse stage and the protostellar stage is illustrated in 
Fig. 1 which compares the radial intensity profiles of the 
pre-stellar core L1689B and the candidate Class 0 object 
L1527. 

We also note that recent large-scale, high-angular resolu¬ 
tion imaging of the p Ophiuchi cloud with the mid-infrared 
camera ISOCAM on board the ISO satellite (Abergel et 
al. 1996) suggests dense cores in clusters are often char¬ 
acterized by very sharp edges (i.e., steeper than p oc 
or p ex r“^), possibly produced by external pressure. (The 
p Oph dense cores are seen as deep absorption structures 
by ISOCAM.) 

These results on the structure of dense cores set con¬ 
straints on the initial conditions for gravitational collapse. 
The starless cores detected in the submillimeter contin¬ 
uum are thought to be descendants of cores without de¬ 
tectable submillimeter emission (which have lower den¬ 
sities) and precursors of cores with embedded low-mass 
protostars. Based on the numbers of sources observed in 
the three groups of cores, the lifetime of starless cores with 
submillimeter emission was estimated to be ~ 10® yr by 
WSHA. So far a dozen pre-stellar cores have been studied 
in detail in the submillimeter, and an inner flattening of 
the radial density profile has been found for all of them. 
Therefore, the statistics indicate that the transition from 
flat to steep inner density profile occurs on a timescale 
shorter than ~ 10® yr, consistent with a phase of fast, pos¬ 
sibly supersonic collapse (Ward-Thompson et al. 1997). 
Ambipolar diffusion models of core formation do predict 
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Fig. 1. Azimuthally averaged flux density profiles of the 
pre-stellar core L1689B (a) and of the candidate Class 0 pro¬ 
tostar L1527 (b) (adapted from AWM96 and Motte et al. 
1996, respectively). Model profiles that fit the data are shown 
as dotted lines. The L1689B model is isothermal and has 
p(r) (X for r < 4000 AU (i.e., 6 < 25") and p(r) oc r~^ 
for r > 4000 AU; the L1527 model has T{r) oc and 

p(r) oc These model profiles result from a full simulation 

of the continuum dual-beam mapping technique (see Motte et 
al. 1996), including a convolution with the beam of the tele¬ 
scope which is shown as dashed lines. Note how flat the L1689B 
profile is compared with the LI527 profile for 6 < 25". The 
apparent flattening of the L1527 profile for 9 f, 5" is entirely 
attributable to the finite resolution of the telescope. 


that the time spent from the critical to the singular state 
is relatively short (typically 10® yr in the models of 
Mouschovias and co-workers). Thus, in itself, the fact that 
most starless cores are observed in the more long-lived 
state when the central density profile is fairly flat does not 
contradict these models. However, close comparison with 
the parameter study presented by Basu & Mouschovias 
(1995) shows that the above timescale estimates agree 
roughly with the magnetically-supported core models only 
if the initial cloud is highly subcritical (Ward-Thompson 
et al. 1997), which implies rather high values of the (static) 
core magnetic field (~ 80 pG in the case of L1689B - see 
AWM96). Upper limits to the field in cloud cores obtained 
through Zeeman observations are typically lower than that 
(e.g. Crutcher et al. 1996, Troland et al. 1996), which sug¬ 
gests the observed dense cores cannot all be supported by 
a static magnetic field and cannot all have formed by am¬ 
bipolar diffusion alone. The same conclusion is indepen¬ 
dently supported by a statistical analysis of the observed 
core shapes which indicates that the overall morphology 
is prolate for most cores (Myers et al. 1991; Ryden 1996). 
This tends to go against ambipolar diffusion models which 
produce oblate cores, flattened in the direction of the mean 
magnetic field. 
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Therefore, although a more complete study of the struc¬ 
ture of starless cores would be desirable, we take the 
present observational evidence to suggest that, in some 
cases at least, the initial conditions for rapid protostellar 
collapse are characterized by a centrally flattened density 
profile and differ significantly from a singular isothermal 
sphere. 


2.2. Decline of outflow power with time 

Class 0 protostars tend to drive highly collimated or “jet¬ 
like” CO molecular outflows (see review by Bachiller 
1996). The mechanical luminosities of these outflows are 
often of the same order as the bolometric luminosities 
of the central sources (e.g., AWB93). In contrast, while 
there is growing evidence that some outflow activity exists 
throughout the embedded phase, the CO outflows from 
Class I sources tend to be poorly collimated and much 
less powerful than those from Class 0 sources. 

In an eflbrt to quantify this evolution of molecular out¬ 
flows during the protostellar phase, Bontemps et al. (1996 
- BATC) have recently obtained and analyzed a homo¬ 
geneous set of CO(2-l) data around a large sample of 
low-luminosity {L^oi < 50 A©), nearby {d < 450 pc) em¬ 
bedded YSOs, including 36 Class I sources and 9 Class 0 
sources. The results show that essentially all embedded 
YSOs have some degree of outflow activity, suggesting the 
outflow phase and the infall/accretion phase coincide. This 
is consistent with the idea that accretion cannot proceed 
without ejection and that outflows are directly powered 
by accretion (e.g., Ferreira & Pelletier 1995). 

In addition, Class 0 objects are found to lie an order of 
magnitude above the well-known correlation between out¬ 
flow momentum flux and bolometric luminosity holding 
for Class I sources. This confirms that Class 0 objects dif¬ 
fer qualitatively from Class I sources, independently of in¬ 
clination effects. On the other hand, as can be seen in 
Fig. 2, outflow momentum flux is roughly proportional to 
circumstellar envelope mass in the entire BATC sample 
(i.e., including both Class I and Class 0 sources). Bon¬ 
temps et al. argue that this new correlation is independent 
of the Aco-Aboi correlation and most likely results from a 
progressive decrease of outflow power with time during the 
accretion phase. 

Since many theoretical models of bipolar outflows (e.g. 
Shu et al. 1994, Ferreira & Pelletier 1995, Fiege & Hen¬ 
riksen 1996a,b) predict a direct proportionality between 
accretion and ejection, BATC further suggest that the ob¬ 
served decline in outflow power reflects a corresponding 
decrease in the mass accretion/infall rate: In this view, 
Afacc would decline from ~ 10 “® M 0 yr“^ for the youngest 
Class 0 protostars to ~ 10“^ M 0 yr“^ for the most evolved 
Class I sources and most active T Tauri stars. 
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Fig. 2. Outflow momentum flux Fco versus circumstellar en¬ 
velope mass Menv for a sample of nearby Class I (filled circles) 
and Class 0 (open circles) YSOs (Bontemps et al. 1996). (The 
arrows on the data points indicate upper limits.) The ‘best 
fit’ Fco-Menv correlation is plotted as a dashed line. This dia¬ 
gram shows a clear decline of outflow power from Class Os to 
Class Is. 


3. Relevant theory 

We propose that the decrease of accretion rate seen by 
BATC during the protostellar phase (see Sect. 2.2 above) 
is a direct consequence of the non-singular initial condi¬ 
tions indicated by the pre-stellar core observations sum¬ 
marized in Sect. 2.1. Our fundamental thesis in this article 
is that supersonic core formation and a subsequent core 
accretion flow may be adequately described as a zero pres- 
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sure flow that evolves from a more general density profile 
than the r~^ profile at a pre-core formation epoch (t < 0). 
We use analytical methods to elucidate the physical ori¬ 
gin of a sharp decrease of M^co during the early accretion 
phase (i.e., just after the completion of stellar core for¬ 
mation). We show that the key to such a decrease is the 
existence of a transition between a fiat inner region and 
an outer inverse square ‘envelope’ in the initial density 
profile: the inner region collapses first to form a hydro¬ 
static stellar core, and the remaining cloud material then 
accretes supersonically onto a non-zero point mass. 

3.1. Overview of known spherical collapse solutions 

Here, we assemble some relevant non-linear studies of 
isothermal protostellar collapse in spherical geometry. The 
isothermality seems justified physically at early (e.g., pre- 
stellar) stages (note this does not imply the dynamical 
importance of pressure), but the assumption of spherical 
symmetry implies the neglect of rotation, magnetic fields 
and dynamical turbulence. 

Mathematically, the paper by Whitworth & Summers 
(1985 - WS85) concludes the series of original studies on 
isothermal, self-similar accretion flows by Larson (1969), 
Pension (1969), Shu77 and Hunter (1977). However physi¬ 
cally the debate continues, based essentially on the failure 
of either the Larson-Penston solution or the Shu solution 
to describe the accretion rate as found in numerical simu¬ 
lations (e.g. Hunter 1977, FC93) after the core-formation 
epoch. The Shu solution describes the asymptotic be¬ 
haviour satisfactorily, while the Larson-Penston solution 
seems more relevant near core formation although it is 
not perfectly so. 

The major contribution of WS85 was to show that 
a parametrically stable infinity of ‘acceptable’ similarity 
solutionaj exist rather than merely the previously known 
Larson-Penston (L/P) infall and the Shu (S) implosion 
behind an outward going expansion wave. ‘Acceptable’ in 
this context implies that the solutions have reasonable ini¬ 
tial conditions inside an inward going compression wave 
in general (that may be initiated by an external distur¬ 
bance) and describe supersonic infall after the passage of 
this wave. True free-fall accretion is established only after 
this wave has reached the centre of the system at t = 0. 
It develops behind the rebounding expansion wave. 

WS85 classified all possible isothermal similarity so¬ 
lutions in a bounded two-dimensional space. Briefly, 
this continuum of solutions spans a wide range from 

^ Similarity solutions are functions defined by ordinary differ¬ 
ential equations which nevertheless represent solutions to the 
full partial differential system of the physical problem. They 
also have the advantage of removing initial and boundary con¬ 
ditions to infinity. The physical justification for their use is 
the empirical fact that they arise as ‘intermediate asymptotes’ 
in various regimes of more complicated flows (Barenblatt & 
Zefdovich 1982). 


‘gravity dominated’ solutions (corresponding to the 
‘band 0/band 1’ of WS85) to ‘sound-wave’ dominated so¬ 
lutions (corresponding to high ‘band’ numbers in WS85). 
The sound-wave (or pressure) dominated solutions are 
those usually studied, which depend on the classical sim¬ 
ilarity variable x = r/at (e.g., Shu77). The S solution is 
itself a limiting sound-wave dominated case (‘band oo’ in 
WS85) that develops from an initial high density singular 
isothermal equilibrium. 

One might think that the gravity dominated solutions 
should be determined in general by the form of the den¬ 
sity or mass profile at some epoch in the region of su¬ 
personic flow, rather than by the isothermality as in the 
above family. We propose a more general type of similar¬ 
ity (see Appendix A and below) that arises from a given 
density profile in a region of supersonic flow. This idea 
was also hinted at in Hunter (1977) in the course of his 
explanation of the absence of shocks in the supersonic re¬ 
gion prior to the arrival of the compression wave at r = 0. 
Hunter observes that gravity provides a non-propagating 
local acceleration by ‘action at a distance’. Formally, the 
flow ceases to be hyperbolic in its causal structure, be¬ 
coming elliptic instead. 

Such a dynamical change in self-similar symmetry was 
considered by Henriksen (1994 - hereafter H94) to argue, 
on the basis of pressure-free, analytical calculations that 
non-constant accretion rates might be expected. We shall 
exploit this idea in the present paper after an assessment 
of the pressure dominated self-similar solutions. 

It is difficult to isolate precisely the epoch at which the 
change in self-similarity type occurs. We know from the 
simulations by Hunter (1977) and recently by FC93 that 
classical self-similarity arises in the course of the collapse. 
The flow is generally found to agree best with the L/P 
type solution, but the agreement is not perfect, since the 
accretion rate found by these numerical studies is not the 
constant value to be expected if the self-similar solutions 
persist at the supersonic core formation epoch {t = 0). 
In fact the accretion rate passes through a peak and then 
settles only gradually (if the cloud is sufficiently large) 
to the Shu value (see Fig. 3 of FC93). The nature of this 
transitory peak was left somewhat unclear by FC93 due 
to their problems of numerical resolution at t = 0, but we 
believe its nature to be vital to the understanding of the 
outflow observations summarized in Sect. 2.2. 

3.2. Justification of our pressure-free analytical approach 

3.2.1. Neglect of magnetic field and rotation 

We neglect both magnetic field and rotation in our 
discussion. Mouschovias and co-workers (Mouschovias 
1995 and references therein) have demonstrated that a 
magnetically-supported, subcritical cloud evolves natu¬ 
rally to a centrally peaked, flattened density profile during 
a quasi-static phase of ambipolar diffusion. This mecha- 
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nism can successfully account for the formation and struc¬ 
ture of some isolated pre-stellar cores. The quasi-static 
phase ends with the formation of a magnetically supercrit¬ 
ical core that collapses in a dynamical but magnetically- 
controlled fashion, while the bulk of the cloud remains 
magnetically supported. In Mouschovias’ models, the ‘dy¬ 
namical’ contraction of the supercritical core is almost al¬ 
ways subsonic and evolves only gradually towards free fall, 
which is not attained up to densities as large as ~ 10®- 
10^® cm“^ (see Basu & Mouschovias 1995 and Fig. 2d of 
Ciolek & Mouschovias 1994). This corresponds to a rather 
gentle process of hydrostatic stellar core formation taking 
significantly longer than a free-fall time, i.e., typically on 
the order of 10® yr in published models from the time the 
parent pre-stellar core becomes magnetically supercritical 
at a central density ~ 10® cm“®. 

However, as pointed out in Sect. 1 and Sect. 2.1, some 
cloud cores are apparently not magnetically supported 
and appear to evolve faster than predicted by ambipolar 
diffusion models. Furthermore, even in the case of initial 
magnetic support, superalfvenic implosion of cloud frag¬ 
ments is possible, for instance induced by a shock (e.g. 
Mouschovias 1989, Tomisaka, Ikeuchi, & Nakamura 1988, 
Tomisaka 1996). This may well be a dominant process 
in regions of multiple star formation. Thus, we believe 
the collapse can become vigorous and supersonic (super¬ 
alfvenic as the case may be) relatively early on, at least 
in a small central mass (see the next section). This is the 
hypothesis we will make in the remainder of the paper. 

Rotation and magnetic fields are thought to be essen¬ 
tial to the launching of the bipolar outflows that are invari¬ 
ably found to be associated with protostars (e.g. Konigl 
& Ruden 1993, Henriksen & Valls-Gabaud 1994). This ar¬ 
gues in favour of their importance, at least when the outer 
parts of the initial cloud core fall in. It is likely that they 
are also important to the formation and evolution of the 
flattened inner parts. 

Thus we cannot expect to include the whole picture 
self-consistently in our subsequent arguments which ignore 
magnetic fields. However if gravity is the main player in 
the core after a ‘point of no return’ is passed, then our 
present approach remains crudely justified, if incomplete. 

3.2.2. Lagrangian treatment of Isothermal Self-Similarity 

As discussed by Henriksen (1989 - H89), spherically sym¬ 
metric, isothermal, self-similar collapse can be very con¬ 
veniently treated when one uses a Lagrangian formal¬ 
ism. This is recalled briefly in Appendix A for the sim¬ 
plest cases. Here, we merely want to show that such a 
formalism provides continuity with the supersonic treat¬ 
ment to be found in Sect. 3.3. The advantage of the 
Lagrangian formulation is that self-similar solutions are 
fully described by the generalized ‘Friedmann, Lemaitre, 
Robertson-Walker’ (FLRW) differential equation (Eq. A2 
in the limit of Eqs. A4 and A5): 


(^) ( 1 ) 

which expresses the conservation of energy. In this equa¬ 
tion, ^ = A/2GM(r)/r® {t — to) is the self-similar vari¬ 
able, where M{r) is the mass distribution (i.e., the mass 
contained within the sphere of radius r as a function 
of r) at some reference epoch t = to (the reason for 
this time offset will become apparent later). The func¬ 
tion S{^) = R{t,r)/r gives the normalized position at 
time t of the shell that was initially {t = to) at r. Fi¬ 
nally, a® = air/GM{r) where as is the isothermal sound 
speed, H = \np/po where Poir) is the density profile at 
t = to, and k = E{r)r/GM{r) = ±1 or 0 is the normal¬ 
ized total energy {E{r) being the specific total energy of 
the shell initially at r). The solutions describing a collapse 
from rest have k = —1, and those with zero total energy 
have fc = 0 (we omit other cases for simplicity). 

The present formulation is unique in that it allows for 
a transition to a purely self-gravitating flow wherein the 
pressure is negligible (formally as 0) before the ingoing 
compression wave establishes the ‘attractor’ at t = 
0 [Blottiau, Bouquet, & Chieze (1988) were the first to 
use this terminology]. In order that a® ^ 0 at t < 0 
in an isothermal flow, it is clear that M(r)/r must tend 
to a large value as r ^ r^, where rjv is the Lagrangian 
radius of this core or nuclear region. This can happen if 
a sufficiently dense mass with a flattened density profile 
forms at the centre of the system at t < 0. This region 
should eventually collapse nearly homologously to form 
a central point mass so long as the inflow is supersonic 
(or super fast-Alfvenic as the case may be) and pressure 
support is negligible. 

Assuming that some fraction of the core does collapse 
in the free-fall time tff (N) , we begin our zero pressure cal¬ 
culations at the onset of this rapid collapse. We define this 
starting point as being to = —tff{N), for consistency with 
the standard pressure wave connected solutions. (With 
this definition, stellar core formation is complete at t = 0, 
which also marks the beginning of the accretion phase.) 
When the parameter a® ^ 0, one has made the transition 
of course from an hyperbolic system to an elliptic one. 
In other words, there is a wave-mitigated causality when 
pressure is important, and action at a distance (in New¬ 
tonian gravity) in the zero pressure limit (see also Hunter 
1977). This is discussed briefly in Appendix A. 

3.2.3. The case for starting from strictly flat inner density 
profiles 

The density profiles observed in pre-stellar clouds 
(Sect. 2.1) are not as centrally flat as are some of the 
self-similar solutions of WS85, nor indeed of the magneto¬ 
static cores of Basu & Mouschovias (1995). This is not 
very surprising since a flat profile indicates a relatively low 
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pressure gradient so that the hydrostatic equilibrium re¬ 
quired for long lived clouds is not possible. In this state the 
cloud is in fact approaching the super-critical state. The 
real question then is whether the collapse proceeds in gen¬ 
eral by some part of the material passing through a phase 
of nearly zero pressure gradient and collapsing homolo¬ 
gously, as the band 0 sound wave dominated self-similar 
solutions suggest (see Figs. 5, 6, 7 of WS85). Interestingly 
enough, the numerical calculations by FC93 that begin 
from non self-similar conditions show a velocity distribu¬ 
tion in the critical Bonnor-Ebert sphere at early times that 
would tend to flatten the density profile (see e.g. the dis¬ 
cussion in the second paragraph of their Sect. 2.3) while 
mildly sub-critical cases show oscillations as are found 
in the ‘bound’ (band 1 and higher) self-similar solutions. 
FC93 interpret their maximum velocity peak in terms of 
a flat interior density distribution that passes abruptly to 
a r~^ profile. In fact, it is easy to show by linearizing the 
equations that an initially flattened density profile will 
first tend to evolve towards a still more flattened config¬ 
uration on a free-fall timescale set by the cloud central 
density. Physically this is due once again to the velocity 
profile of the fastest growing linear mode that is consis¬ 
tent with zero velocity at the cloud edge and at the cloud 
centre, since this is such that the outer layers overtake the 
inner ones. The linear calculation can only be suggestive, 
but in the generalized solutions of Blottiau et al. (1988) 
as well as in the band 0 solutions of WS85, flattened cen¬ 
tres and shells (which would also collapse to give a point 
mass) appear. Moreover during the ‘quasi-static’ evolution 
followed by Basu & Mouschovias (1995), flattened super¬ 
critical cores appear. That a portion of this core collapses 
homologously is however less clear from their calculations 
(depending on the ion-neutral collision time). 


3.3. Accretion phase history 

For the reasons given in Sect. 3.2 above, because of the 
nature of the abundant non-linear self-similar solutions 
(WS85), and because of the suggestive numerical study of 
FC93, we now proceed to calculate an accretion history 
for a self-similar cloud based on the following hypothesis 
of pre-collapse evolution. 

In agreement with the observations, we suppose that the 
pre-collapse cloud core evolves quasi-statically, and devel¬ 
ops a dense central region with a nearly flat density pro¬ 
file, until an external disturbance initiates the collapse. 
We further suppose that, after a transitional period (lead¬ 
ing to, e.g., magnetic decoupling), the flat region is free 
to collapse suddenly, essentially homologously due to the 
absence of a pressure gradient. It is possible that further 
flattening of the profile occurs during the collapse. In any 
case the beginning of this free collapse is taken once again 
to he at t = to = —tff{N), so that the completion of the 
core collapse is at t = 0, where tff{N) is the free-fall time 
of a uniform sphere from rest, defined explicitly below. 


For all times t > to we shall assume that the evolution is 
of the gravity-dominated self-similar form, given uniquely 
by the density profile at the t = to epoch (see Appendix 
A). This is all that is required for the purposes of our cal¬ 
culation, but the flow has probably ceased to be pressure 
dominated even before t = to. 

An important point to realize is that the accretion rate 
at t = 0 will not yet be that appropriate to the exterior 
of the cloud, since this will be for some time dominated 
by the collapsed core mass rather than by the self-gravity 
appropriate to the outer density profile. The subsequent 
transition to a constant asymptotic accretion rate hap¬ 
pens as sufficient mass from the outer cloud is accreted to 
dominate the initial core mass. 


3.3.1. General Formulation 


The equations include Eq. (1) above, plus those of Ap¬ 
pendix A and of H89, H94 but with = 0. We carry 
out the calculations for a family of initial density pro¬ 
files which are all strictly flat for r < rN (see Fig. 3 
below). The density is measured in units of the density 
Pn of the flattened nucleus of the pre-collapse cloud (i.e., 
p = p/pn), the mass is in units of the mass nucleus 
Mat (i.e., M = M/M^) and the radius is in units of 
Tat (i.e., f = r/r^), where tat is the radius of the nu¬ 
clear region just before free-fall collapse. For convenience, 
time is measured from the beginning of the free collapse 
{t = to = —tff{N j) in units of the nuclear region free-fall 
time tff{N) = (37r/32G/9Ar)^^^, that is we introduce the 
dimensionless variable 


_ t-to 

iff {NY 


( 2 ) 


The stellar core is fully formed at r = 1 and consists 
of all shells with Lagrangian labels f < 1. We take this 
epoch to correspond to ^ = 7r/2 since this is appropriate 
to the collapse of a uniform sphere from rest (e.g. Hunter 
1962). In this description, t = 1 also marks the beginning 
of the YSO accretion phase, while dimensionless times in 
the range 0 < r < 1 correspond to ‘isothermal protostars’ 
(i.e., collapsing cloud fragments with no central hydro¬ 
static core - see Sect. 1). 

The introduction of the self-similar variable C (0 ^ C ^ 
7r/2) is partly justified by the fact that the free-fall time 
of the shell with Lagrangian label r is related to the mass 
function M{r) through 


iffY") = 


r"’ . ,,,, / 7 ’‘/ 2 , 


M„) «//(«) a 


Thus, for f > 1 we have ^ < 7r/2 and 


7r/2 



( 3 ) 
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(In the post nuclear formation epoch, we have ^ = 7r/2 
for all f < 1 since these shells have already collapsed in 
the singular nucleus.) 

We shall be interested in those shells which are ac¬ 
creting near the center during the post nuclear formation 
phase and thus which are about to enter the singularity. 
As is well known in cosmology (FLRW models, see H89), 
for such shells a collapse in the absence of significant pres¬ 
sure gradients can be described by Eq. (1) with A: = 0 (and 
ctg =0), for which the solution is 

/o \ 2/3 

, ( 4 ) 


It is perhaps worth repeating this equation in the form 
that it assumes with restored dimensions, namely 


M*(t) 


7r2 r3(t) 
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(70 


Conversely, the time at which the shell initially at f 
accretes onto the central singularity is 



As for the total mass of the accreting circumstellar enve¬ 
lope at time r, it is simply Menvir) = Mcloud - A7*(r). 


for ^ < 7r/2. 

We note immediately that t = I when ^ = 7r/2 for 
the outer shell of the nucleus (f = 1), which is the instant 
when the nucleus has fully formed. 

To obtain the relationship between the real position of 
a spherical shell R{t, r) and its Lagrangian or comoving 
label r we must use R = rS{^), which with Eq. (4) gives 


X 

7r/2 



( 5 ) 


We see that R = 0 when ^ = 7r/2 as expected. Note how¬ 
ever that nevertheless the shell label r will be greater than 
or equal to I in our calculations, since we shall be study¬ 
ing spheres whose initial radii are greater than or equal to 
tn- 

Our goal is to calculate the accretion rate history and 
the mass remaining in the envelope surrounding the col¬ 
lapsed nucleus, since both of these quantities are suscepti¬ 
ble to observation. The accretion rate is given at each 
spherical radius R and time t by = AttR"^ pdtR. If we 
use our various definitions and measure distance and den¬ 
sity in units as above, and velocity in units of the nuclear 
free-fall velocity vn = yj2 GMN/tn we obtain 


M = _ — _ 

acc — A 9 

iirrf^PNVN 




dS M{f) 


( 6 ) 


where S is given by Eq. (4) and p follows from Eq. 
in the form 


P = 


Po{r) 


^3(1 


din5 \ ’ 

d In^ d\nr ^ 


(A3) 


where po is the density profile at t = to- 

Since we are mainly interested in calculating the accre¬ 
tion rate near the centre of the post-collapse configuration 
where R/r —> 0, we may let ^ tend towards ^ = 7r/2 in 
Eq. (6). 

The mass M* of the central collapsed object accreting 
the shell labelled f at the time r is also easily obtained by 
setting ^ = 7r/2 in Eq. (3) : 


M*(t) 


r'3(r) 




( 7 ) 


3.3.2. Two-component Power-Law Model 

To obtain explicit expressions for and Menu , we must 
assume an explicit initial density profile. We use a multi- 
component profile of the form (in the units introduced in 
Sect. 3.3.1) 

p'o = 1 < r < fb, (9) 

= 7>r'b, (9') 

n 

and Po = 1 for f < 1, where r?, is a characteris¬ 
tic boundary radius, Di is a positive number character¬ 
izing the power-law slope of the cloud ‘envelope’, and 
0 < II 2 < 2/3 characterizes the radial profile of the ‘edge’ 
of the model pre-stellar core (see Fig. 3). This form gives 
a weak discontinuity in the density profile between the 
flat nucleus and the envelope (the density remains 

continuous but its derivative is discontinuous at r = r^r). 
This discontinuity breaks the density self-similarity but 
not the dynamical self-similarity as discussed in H89. In a 
pressure-free epoch there is no reason why such a discon¬ 
tinuity might not be preserved, if ever it were to be estab¬ 
lished by an inward propagating wave for example. Such 
weak discontinuities often arise during the adjustment of 
a gas to a change in the boundary conditions (e.g. Landau 
& Lifshitz 1987) and this may well be the manner in which 
the collapse is initiated. In this subsection, we will keep Di 
and D 2 as free parameters in order to obtain general for¬ 
mulae and model a wide range of accretion histories, but in 
practice we will use £>1 = 1 (r“^ envelope) and £>2 =0.1 
when comparing with observations (Sect. 4). The rationale 
behind such a two-component, idealized density profile is 
that it approximates the observed structure of pre-stellar 
cores reasonably well (cf. WSHA, AWM96, Abergel et al. 
1996, and Sect. 2.1) and can account for the necessarily 
finite radius of influence of protostars forming in clusters 
(since £>2 < 2/3 the model pre-stellar core has a steep 
edge and a finite total mass, see below). With this den¬ 
sity profile, the mass function M(r) that appears in the 
self-similar variable ^ (see Appendix A and H89) in the 
post ‘nuclear collapse’ domain (^ > 0) is podr -\- Mm, 
where the nuclear mass is Mm = (47r/3)/9Arr^. Thus, the 
mass distribution in units oi Mm takes the simple form 
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Fig. 3. Density vs. radius (in log-log format) for an idealized 
pre-stellar cloud with a flat nucleus surrounded by a power-law 
envelope and a steeper power-law outer boundary 


M(r) = - 2)/{3Di - 2), 1 < f < h, 

( 10 ) 

QTO — 2/t3l) ~ 

M(f) = M,+ - 1], f > n, 

6 D 2 — 2 Tb 

( 10 ') 

where Mb = 2 /_Di) _ 2 )/( 3 T)^ _ 2 ) is the mass 

enclosed within the boundary radius fb- Since D 2 < 2/3, 
the mass function converges at large f towards Mdoud = 
Mb + V(2 — 3112 )] • (Note that the case D 2 = 

2/3, which we do not include explicitly here, is logarith¬ 
mic.) 

In this case, the time function r(f) giving the time of 
arrival in the central hydrostatic core of the shell initially 
at f (see Eq. 8) can be easily inverted and has simple 
power-law asymptotic forms: f(r) ~ ^ 

before accretion of the shell labeled fb, and 
r(r) ~ ^cioud ^ later times. Figure 4 displays this 

r(r) relationship for three initial density profiles which all 
have a sharp outer boundary (fb = 10 and D 2 = 0.1). The 
special = 1 case (r“^ ‘envelope’, which we will use in 
practice) is particularly instructive since it shows that a 
singular pre-stellar core would collapse according to 
the linear relation f = \/3t or r(t) = vn (t — to) in 
restored dimensions. This makes possible a direct compar¬ 
ison with Shu77’s expansion wave solution which predicts 
r(t) = Rinfl2 = (as/2)t (here denotes the head of 
the expansion wave at time t; note that to = 0 in Shu’s 
solution). In a sense, Fig. 4 illustrates that, during the 
accretion phase, the idealized solutions explored in this 


Fig. 4. Radius of the shell accreting onto the central singular¬ 
ity as a function of time for three values of Di corresponding 
to r~^ (dotted curve), r~^ (solid curve), and (dash-dotted 
curve) envelopes respectively in the initial density profile. A 
sharp outer edge (fb = 10 and D 2 = 0.1) has been assumed 
in all three cases. The dashed line shows the linear relation 
f = expected for a singular r~^ sphere. By dehnition, 

T — 1 corresponds to the time at which the outer shell of the 
flat nucleus reaches the central singularity 


paper are characterized by an ‘inside-out’ collapse behind 
an expansion wave reminiscent of the Shu solution. 

Substituting Eqs. (10) and (10’) into Eq (6), we obtain 
the desired accretion rate history at ^ = 7r/2 explicitly as 


M 


acc 


2 r(s-2/D,) 3j4^f(3-2/ttD _ 2 

37r 7- f(3-2/r>i) _ I 


( 11 ) 


as long as the expansion wave is still within fb- At 
later times, the accretion history may be obtained from 
the relation: 


M 


acc 


4 f^^^f^^ 2/n2) M(f) 


(IT) 

where D is defined by 1/D = (I/D 2 ) — (l/Di). Based on 
these formulae. Fig. 5 displays the predicted accretion rate 
as a function of time for the same three model pre-stellar 
profiles as in Fig. 4. The corresponding relationships be¬ 


tween Mace and Mgnv are shown in Fig. 6. 

Equation (11) implies that the accretion rate diverges 
at f = 1 (i.e., r = 1) and asymptotically approaches a 


power-law dependence at large f (i.e., large r), Mace oc 
t(3^i- 3) if Di > 2/3, and Mace <X rl^D^-4/3D^) ^ 
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TIME, T (in units of Tff^) 


Fig. 5. Accretion rate as a function of time for three values of 
_Di corresponding to r~^ (dotted cnrve), r~^ (solid cnrve), and 
(dash-dotted curve) envelopes respectively in the initial 
density profile. A sharp outer edge (f(, = 10 and D 2 = 0.1) has 
been assnmed in all three cases. The dashed horizontal line 
gives the asymptotic value of the accretion rate obtained in 
the case an r~^ initial density profile 


2/3 (see Fig. 5). (Here again, the case Di =2/3 is log¬ 
arithmic.) In particular, Mace tends to the asymptotic 
value —•\/3.(2/7r) at large r when Di = 1 (r“^ envelope) 
(see also Appendix B). As was observed by H94 and FC93, 
this asymptotic limit is only achieved if the envelope is 
snfhciently large (i.e., typically fb ~ 4). Otherwise the ac¬ 
cretion rate associated with stellar formation is truncated 
abruptly as the surface of the protostellar cloud falls in. 
Terminal accretion may continue after this stage as noted 
in H94 but with a character that depends on the density 
profile in the environment. We try to simulate this late 
effect here with a steep D 2 = 0.1 edge (see Eq. 9’). 

The accretion peak at r = 1 is due to the close approach 
of the density profile to being flat near f = 1. Physically, 
if the density profile is strictly flat, two distinct spher¬ 
ical shells fall to the centre simultaneously, so that the 
frequency of shell arrival (which is the accretion rate) is 
infinite. It is in fact easy to show from the preceding the¬ 
ory that the criterion for a shell at r and another at r + 6r 
to be accreted at i? = 0 at a transitory peak rate is indeed 
that 


d\nM{r) 

dlnr 


( 12 ) 


In fact, if the density profile is not perfectly flat at some 
point Tat, e.g., Po oc locally where e is small and pos¬ 
itive, then the logarithmic derivative of Eq. (12) attains 
the value 3 — 2e, and one finds that the accretion rate for 
these shells will have a peak amplitude (in the usual units 
at rjv) 

peak O 

Mace = -• ( 13 ) 

Tre 



Fig. 6. Absolute value of the accretion rate as a function of 
envelope mass for the same three initial density profiles as in 
Fig. 4 and Fig. 5. Here again, the dashed horizontal line gives 
the asymptotic value of the accretion rate predicted in the case 
an initial density profile 


4. Discussion: Implications for protostellar evolu¬ 
tion 

What we should retain from the previous analytical cal¬ 
culations is the association of centrally peaked, flattened, 
pre-stellar density profiles with a transitory, vigorous, 
post nuclear formation accretion epoch. We here suggest 
that this energetic accretion phase coincides with Class 0 
protostars, which would explain their unusually power¬ 
ful jets compared to the more evolved Class I YSOs (see 
Sect. 2.2). 

4.I. Predicted versus observed accretion histories 

We proceed to compare the model predictions obtained 
in Sect. 3.3 above with the Fco c/Tboi versus Me„v/Fboi°'® 
diagram corresponding to the entire BATC sample. This 
diagram, which is almost free of any luminosity or distance 
effect, should mainly reflect the evolution of the accretion 
rate Macc- Indeed, as discussed by AWB93 and BATC, 
the M_/Lboi° ® ratio can be used as a practical, quan¬ 
titative indicator of embedded YSO evolution, which, to 
first order, should track the ratio of envelope to 
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stellar mass and thus decrease with protostellar age.0 As 
to the CO outflow momentum flux Fqo, it is related to 
the accretion rate M^cc by 

FcO = [/e„t (Mw/Macc) K] X M_, (14) 

where /e„t is the entrainment efficiency, is the 
mass-loss rate of the underlying driving wind or jet, 
and V„ is the wind velocity. As discussed in BATC, 
/e„t (M„/Maoc) K. is unlikely to vary much during the pro¬ 
tostellar phase, so that Fco should reflect the variations 
of the accretion rate However, because as a gen¬ 

eral trend more massive cores tend to form more massive 
and more luminous YSOs associated with more energetic 
outflows (e.g. Saraceno et al. 1996), F^o also depends on 
luminosity or stellar mass. For this reason, we use the 
quantity F^o instead which should be less sensitive 

to the initial cloud mass/density. To the extent that there 
is no significant luminosity evolution for a given proto¬ 
stellar source (see Kenyon & Hartmann 1995 and BATC), 
Fco c/Fboi should more homogeneously trace the intrinsic 
temporal variations of M^cc than Fco- 

The Fco c/Fb„i versus diagram observed 

by BATC, which is shown in Fig. 7, suggests that the ac¬ 
cretion rate of Class 0 sources is a factor of ~ 10 larger on 
average than that of Class I sources (see also Sect. 4.3.1 
of BATC). By comparison with Fig. 6, it is tempting to 
identify the Class 0 phase with the short period of vig¬ 
orous accretion predicted by the theory of Sect. 3.3 after 
the formation of the hydrostatic stellar nucleus, and the 
Class I phase with the longer period of constant M„c pre¬ 
dicted at late times when Fi = 1 (r“^ envelope). 

To make quantitative comparisons, we use an in¬ 
dependent observational constraint on the average val¬ 
ues of the accretion rate during the Class I and Class 0 
phases. The lifetime of Class I sources is measured to 
be At/ ^ 2 X 10® yr in clouds forming mainly low-mass 
stars such as p Ophiuchi and Taurus (e.g. Wilking et al. 
1989, Kenyon et al. 1990, Greene et al. 1994). The life¬ 
time of Class 0 sources is estimated to be an order of 
magnitude shorter, Atg ~ 10"’’ yr (see AM94 and Barsony 
1994). In this paper, we adopt the conceptual dehnitions 
of the Class 0 and Class I stages introduced by AWB93 
(see also AM94). The Class 0 stage corresponds to proto¬ 
stars which have developed a hydrostatic stellar core but 
still have > M*, i.e., objects which have accreted 

less than a half of the total cloud mass; in our model, the 
formation of the central hydrostatic object ends when the 
flat pre-stellar nucleus has fully collapsed, thus Class 0 
sources have M// < M, < Mcioud/2. The Class I stage 
corresponds to embedded sources with smaller envelope 

^ The envelope masses (Me„v) used in the diagrams shown by 
BATC were derived from 1.3 mm maps by integrating the ob¬ 
served dust continuum emission over a region ~ 1' in diameter 
(corresponding to ~ 10"’ AU at d = 160 pc). A mass opacity 
ki .3 = 0.01 cm^g“’ was assumed (cf. AM94). 



Fig. 7. Fcoc/L boi (dimensionless) versus A7env/Tboi (IHenv 
and I/boi in units of Mq and Lq) for the sample of Class 0 
(open circles) and Class I (filled circles) sources studied by 
Bontemps et al. (1996). (The arrows on the data points indicate 
upper limits.) Fco c/Lboi can be taken as an empirical tracer of 
A7acc provided that there is no significant luminosity evolution. 
Me„v/Fboi°'® is an evolutionary indicator which decreases with 
protostellar age. The eye-fitted solid curve shows the accretion 
rate history predicted by the model of Sect. 3.3.2 assuming 
Di = 1, D 2 ~ 0.1, and ft, = 1.6 (see text). The dashed line 
shows the corresponding asymptotic value of Fcoc/Lboi for 
comparison with Fig. 5 and Fig. 6. 

masses: Me„v < M„. We also define the boundary between 
Class I and Class H sources by M„ = 0.99 x Mdoud, i.e., 
Afenv = 0.01 X Mcioud- TMs allows for residual (disk) ac¬ 
cretion during the (Class H) pre-main sequence phase: a 
1 Mq star may typically accrete up to the mass of the 
minimum solar nebula, i.e., 0.01 Mq, at the Class H stage 
(0.01 Mq is indeed the median circumstellar disk mass 
measured around Class H T Tauri stars - Beckwith et al. 
1990, AM94). 

Within the framework of the model introduced in 
Sect. 3.3.2 the mass function M(f) is easy to invert, and 
given the above definitions the lifetime ratio Atj/Ato can 
be calculated analytically as a function of Mdoud or f/: 

Ati/Ato = (r/ - To )/(to - 1), (15) 

where 

To = [r{Mdoud/2)fV{Mdoud/2)^^^ (16) 

and 


TI = [f{0.99Mdoud)fViO-^^Mdoud)^^^. 


(16') 
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This model ratio is plotted in Fig. 8 as a function of 
the boundary radius fb, for the same three values of 
Di as in Figs. 4, 5, 6. It can be seen from Fig. 8 that 


M 

O 




Fig. 8. Predicted ratio of the duration of the Class I phase 
Ati to the duration of the Class 0 phase Ato as a function 
of the pre-stellar boundary radius fb within the model of 
Sect. 3.3.2 (see text for formal definitions of the Class 0 and 
Class I phases). Three representative cases are shown: r~^ en¬ 
velope {Di = 2, dotted curve), r~^ envelope (Hi = 1, solid 
curve), envelope (Hi = 2/3, dash-dotted curve). In all 
three cases, a sharp outer boundary (H 2 = 0.1) has been as¬ 
sumed beyond fb- 


the observational constraint Atj/Ato ~ 10 can be sat¬ 
isfied by the Hi = 1 model (standard case) only if the 
pre-stellar boundary radius is small (ff, ~ 1-6), corre¬ 
sponding to a small value of the normalized cloud mass 
[Mcioud ~ 3). Furthermore, the predicted At//Ato ratio 
decreases rapidly towards an asymptotic value of 1 as f/ 
increases (one already has Atj/Ato ^ 2 for fb = 3.5, i.e., 
Mcioud ~ 9). The Hi = 2 case is qualitatively similar. In 
contrast, for sharper initial density prohles such as r~^ 
(Hi = 2/3), the predicted lifetime ratio At//Ato is ~ 10 
for a wide range of fb values. 

Since the mean stellar mass in the BATC sample is 

< M* > ^ 1 Mq, the average accretion rate during the 
Class I phase can be estimated as 

< >/«< M* > /(2At/) ~ 2.5 X 10-6 Mgyr"!. In 
addition, one has 

< M^^c >0 / < M^^c >i ~ At//Ato X (1 ~ -)• (17) 

^cloud 

Therefore, if Mdoud is sufficiently larger than 2, the 
lifetime constraint At//Ato ~ 10 implies 


< Macc >0 / < Afacc >/ ~ 10 , in reasonable agreement 
with the drop of Fco c/Tboi seen in Fig. 7 between Class Os 
and Class Is. 

To reduce the number of free parameters, we will as¬ 
sume Hi = 1 and H 2 = 0.1 from now on. With these 
values, the model profile accounts for the nearly r~^ ‘enve¬ 
lope’ that is observed in most pre-stellar cores just outside 
the flat central region (WSHA, AWM96). It also features 
a steep outer edge, qualitatively similar to those revealed 
by some recent infrared absorption studies of dense cores 
(Dent et al. 1995, Abergel et al. 1996). 

The boundary radius fb is the only parameter of the model 
that we will vary in order to match the observations. As 
discussed above, the lifetime constraint suggests fb ~ 1 . 6 . 
In Fig. 7, we have superposed as a solid curve the accre¬ 
tion history predicted by the model when f/ = 1.6 (and 
Hi = 1 , H 2 = 0 . 1 ) to the Fcoc/L^c .1 versus 
diagram observed by BATC. To scale the model along the 
x-axis, we adopted Mdoud = 1 Mq, based on the average 
value of Me„v for the Class 0 sources in the BATC sam¬ 
ple. The outflow efficiency factor fent{M„/M^cc)V„ was 
adjusted so as to make the model match the typical value 
Fcoc/Lboi ~ 200 observed among the Class I sources of 
the BATC sample. Given an average bolometric luminos¬ 
ity < Lboi > = 2 Lq (BATC), this corresponds to a typ¬ 
ical momentum flux < F/o > ~ 8 x lO”® M 0 km s-^yr”^ 
for Class I sources. At late times the (Hi = 1) model 
predicts a constant accretion rate M^^^ corresponding to 
the standard Shu value (see Appendix B). If we adopt 
Macc ~ 2 X 10-6M 0 yr-^, which is appropriate for a 
T = 10 K cloud and in rough agreement with the value 
estimated above from the lifetime of Class I sources, then 
we find that the best ‘eye-fit’ model has an efficiency 
/e„t {M„/M^cc) V„ ~ 4km s-^. We note that this is a fac¬ 
tor of 4 lower than the efficiency used in BATC which 
corresponds to the canonical values V„ = 150 km 3 -^, 
= 0.1, and /e„t = 1. A possible explanation for 
this difference is that the entrainment efficiency is actually 
lower than 1, /e„t ~ 0.25. In spite of this discrepancy, the 
overall agreement between the shape of the predicted ac¬ 
cretion history and the observations is quite encouraging 
given the simplicity of the present modeling. In particu¬ 
lar, it should be noted that the range of initial Mdoud 
values spanned by the BATC sample is certainly only 
partly accounted for by dividing Fbo by Lboi and plot¬ 
ting Fco c/Tboi versus M,„„/Lboi° ®- 

A potential problem of the H/ = 1 model shown in 
Fig. 7 is that it apparently leads to very short timescales. 
For fb = 1.6, the predicted Class 0 and Class I lifetimes 
are indeed Affi ~ 0.04f//(A) and Atj ~ 0.37tff{N) 
in units of the nuclear region free-fall time. If we take 
the mean density n //2 ~ 2 x lO^ cm-^ observed in the 
flat central region of L1689B (AWM96) as a representa¬ 
tive example, we obtain tff{N) ~ 7 x 10'^ yr, and thus 
Ato ~ 2500 yr and Atj ~ 2.6 x lO"* yr, which are at least 
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a factor of ~ 4 shorter than the observed values of Ato 
and At I. 

However, in actual fact, the collapse is probably less vio¬ 
lent than indicated by our self-similar, pressure-less calcu¬ 
lations. The more realistic numerical calculations of FC93 
and Tomisaka (1996), which fully account for pressure sup¬ 
port (including magnetic support in the latter case), yield 
qualitatively similar accretion histories but quantitatively 
longer timescales. In particular, the core formation time 
since the initiation of collapse is a factor of ~ 4-5 times 
larger than the free-fall time tff{N) in these calculations. 
If we use this factor to roughly estimate the slowing down 
influence of pressure effects, then we find Ato ~ 10^ yr 
and At I ~ 10^ yr, in better agreement with the observed 
lifetimes. 


4-2. Accretion luminosities of Class 0 and Class I sources 

If, as proposed here, Macc is a factor of ~ 10 larger for 
Class 0 sources than for Class I sources, one may naively 
expect the former to have much higher accretion luminosi¬ 
ties than the latter. This would be in apparent contradic¬ 
tion with observations which suggest that Class Os are not 
significantly over-luminous compared to Class Is. For in¬ 
stance, in the embedded YSO sample of BATC, the ratio 
of the average bolometric luminosities for the two classes 
was < Lboi >0 / < Aboi >/~ 1.6. (Note however that the 
typical luminosity of Class 0 sources in a given cloud is 
relatively uncertain observationally due to the rarity of 
these objects.) 

In reality, several factors contribute to reduce the theo¬ 
retical accretion luminosity (Lace = GM-i,Macc/R*) at the 
Class 0 stage, making it comparable to the Class I accre¬ 
tion luminosity: 

(1) The central stellar mass M* is smaller for Class 0 
sources; 

(2) The stellar radius i?* is likely to be larger if is 

• 1 /3 

higher, since one expects the rough scaling i?* oc 
(Stabler, Shu, & Taam 1980; Stabler 1988); 

(3) The amount of accretion energy dissipated in the wind 
can be expected to be larger for Class Os than for Class Is 
(and could be a significant fraction of the total accretion 
energy). 

Ignoring effect (3) for a moment, it is easy to compute 
a rough theoretical estimate for the ratio 

< >0 / < >/. Indeed, 

< > 0 = ^ 


^ Laicc ^ 0 ^ 


1 


G 


('Mcloud!^ 


A^o < -Ra >0 j mn 
Likewise, < >/= 

1 G 


(18) 


To K* ’ 


^cloud 


^ -^acc ^ 


Ati < >i Ju^ioudl^ 


M^dM,. (18') 


Thus, 

< >0 ^ 1 Atl < R^, >I ^ 4 

< >, ^ 3 Ato < i?* >0 ^ 

In these equations, < i?* >o and < i?* >/ are ap¬ 
propriately time-averaged values of the stellar radius at 
the Class 0 and Class I stages, respectively. If the scaling 
Ri, oc is valid, then one predicts ^ ;$ 1.5 for 

10 , which is consistent with existing observations. 
The magnitude of effect (3) is more difficult to assess since 
it depends on the still poorly understood mechanism of ac¬ 
cretion/ejection. In the case of the X-celerator mechanism, 
the fraction of power extracted by the wind roughly scales 
as the ratio R^,/ Rx where Rx is the disk truncation ra¬ 
dius (e.g. Shu 1995). As the accretion rate decreases and 
the centrifugal disk radius increases, we may expect Rx to 
increase, thereby reducing the fraction of power released 
in the wind. The order of magnitude of this effect may be 
another factor of ~ 2 in Lacc- 

In conclusion, the combined effects of (1), (2), (3) 
above are likely to render the luminosities of Class 0 
sources similar to those of Class I sources. 


4-3. Differences between p Ophiuchi and Taurus 

When comparing model predictions with observations in 
Sect. 4.1 we considered all the YSOs of the BATC sam¬ 
ple together, regardless of their parent clouds. However, if 
we consider the Fco c/Lboi versus Me„v/Tboi°’® diagrams of 
p Ophiuchi (Fig. 9) and Taurus (Fig. 10) separately, clear 
differences become apparent. In the diagram of Fig. 9, 
the Class 0 sources of p Oph (VLA 1623 and IRAS 16293) 
clearly stand out as a distinct group characterized by val¬ 
ues of J^coc/Tboi and ® which are both an or¬ 

der of magnitude larger than the corresponding values 
for Class I sources. In fact, this clear contrast observed 
in p Oph between Class 0 and Class I objects was part 
of the original motivation for introducing the Class 0 as 
a new class of YSOs (see AWB93). On the other hand, 
no such contrast is observed in the Taurus diagram of 
Fig. 10, where there is a much better continuity between 
Class 0 and Class I sources. In other words, the Class 0 
candidates of Taurus may merely correspond to “extreme 
Class I” sources. 

Comparison of Fig. 9 and Fig. 10 suggests that p Oph 
and Taurus YSOs follow different accretion histories. We 
suggest that this directly results from differences in ini¬ 
tial conditions. Indeed, an important feature of the ac¬ 
cretion model proposed in Sect. 3.3.2 (see also FC93) is 
that it predicts a significant drop of the mean accretion 
rate from the Class 0 to the Class I stage only if the ra¬ 
dius rjv of the flat inner region in the initial profile is a 
large fraction of the cloud boundary radius rt,, i.e., if ff, is 
small (see Fig. 8). Thus, the observed difference between 
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Ophiuchus 




Fig. 9. Fdo r./ Tboi versus IWenv/ diagram for tire 
sub-sample of Class 0 (open circles) and Class I (filled circles) 
sources observed by BATC in p Ophiuchi. (The arrow indi¬ 
cates npper limits in both variables.) The eye-fitted solid curve 
shows the accretion rate history predicted under the model of 
Sect. 3.3.2 with ff, = 1.5, D\ = 1, and D 2 = 0.1 


Taurus 



Fig. 10. Fcoc/Lboi versus Me„v/Tboi°'® diagram for the 
sub-sample of candidate Class 0 (open circles) and Class I 
(filled circles) protostars observed by BATC in Taurus. (The 
arrows on the data points indicate upper limits.) The eye-htted 
solid curve shows the accretion rate history predicted by the 
model of Sect. 3.3.2 assuming f(, = 5, Hi = 1, and D 2 = 0.1 


Fig. 9 and Fig. 10 could be accounted for, at least quali¬ 
tatively, by the present models if the boundary radius fb 
were significantly smaller in p Ophiuchi than in Taurus. 

Interestingly enough, such a trend in the values of fb 
seems to be borne out by independent (sub)millimeter 
continuum results on the structure of protostellar en¬ 
velopes and dense cores. In Taurus where YSOs form 
in relative isolation, protostellar envelopes are observed 
to be extended over radii ^ 10^ AU (e.g. Motte et al. 
1996). In contrast, p Ophiuchi is a star-forming cluster 
where the fragmentation size scale is significantly smaller 
than 10'^ AU (e.g. Motte et al. 1997), so that the ra¬ 
dius of the ‘sphere of influence’ of a given YSO must be 
rb < 4000 AU. The density profiles of pre-stellar cores are 
always found to flatten out at radii r < tm ^ 4000 AU 
(see Sect. 2.1 and references therein), with no clear dif¬ 
ference between p Oph and Taurus. These observations 
suggest that, if bounded at all, the Taurus cores have a 
large boundary radius fb > 3, while the p Oph cores have 
a smaller boundary radius r;, ^ 1. 

Models constructed with fb = 1.5 and fb = 5, cor¬ 
responding to Mcioud ~ 2.8 and M^ioud ~ 13.9, have 
been superposed on the diagrams of Fig. 9 and Fig. 10, 
respectively. The models shown in the figures predict 

< Mace >0 / < >/~ 4 in Ophiuchus and 

< Mace >0 / < .^acc >/~ 1-4 m Taucus, in reasonable 

agreement with observations. As in Sect. 4.1, the outflow 
efficiency factor /c„t K, was adjusted in order 

to ‘fit’ the observations (by eye). We adopted standard, 
asymptotic accretion rate values of M°^^ ~ 10“® M0yr“^ 
in Ophiuchus and ~ 2 x 10“® M0yr“^ in Tau¬ 

rus (e.g. Adams, Lada, & Shu 1987), and we used av¬ 
erage bolometric luminosities of < Lboi > = 6 L0 and 

< Lboi >= 0.7 L0 for the Class I sources of the two 

clouds, respectively (see BATC). This procedure yields an 
outflow efficiency /e„t V„ ~ 6 km s“^ in p Oph, 

and /e„t (M„/Macc) K, ~ 1.4 km s~^ in Taurus. (Taken at 
face value, these numbers point to a somewhat more effi¬ 
cient ejection process in p Ophiuchi than in Taurus.) 

In p Oph, the model fit of Fig. 9 suggests that the 
Class 0 sources are in the initial phase of vigorous accre¬ 
tion while the Class I sources may be in the terminal ac¬ 
cretion phase. The latter is not surprising since, due to the 
small fragmentation size scale of this cloud, the expansion 
wave characterizing all collapse solutions at t > 0 (Shu77, 
WS85) will reach the boundary of any given protostellar 
envelope/core in a time tb = rb/ag ^ 10^ yr (assuming 

Os = 0.35km s“^), shorter than the typical Class I life¬ 
time. 

In contrast, all the Taurus sources of Fig. 10 appear to be 
in the ‘asymptotic’ phase during which the accretion rate 
is approximately that predicted by the Shu theory. 

The phase of enhanced accretion is thus apparently ob¬ 
servable in p Oph but not in Taurus. We suggest this is 
related to the fact that core collapse is probably induced 
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by the impact of a weak interstellar shock wave in the first 
case (e.g. Vrba 1977, Loren & Wootten 1986) and self- 
initiated by ambipolar diffusion in the second case (e.g. 
Lizano & Shu 1989, Mouschovias 1991). 

We conclude that the standard theory of Shu and co¬ 
workers is roughly adequate in Taurus but incomplete in p 
Oph. In the latter cloud, the accretion scenario advocated 
here provides a better fit to the data. Further support 
to this view comes from the fact that the only candidate 
‘isothermal protostars’ known to date have precisely been 
found in regions of multiple star formation like p Oph, 
in the form of very dense compact clumps only seen at 
submillimeter wavelengths (Mezger et al. 1992a,b; Chini 
et al. 1997). Indeed, while the initial isothermal collapse 
phase is vanishingly short in the standard Shu picture (see 
Sect. 1), it should last for a significant fraction of the total 
collapse/accretion time if rjq jrh is large and if the scenario 
proposed in this paper is approximately correct (see, e.g.. 
Sect. 3.3.1). 

5. Summary and conclusions 

The main points of our paper may be summarized as fol¬ 
lows: 

1. The radial density gradient of pre-stellar cores is typ¬ 
ically flatter than p{r) oc r~^ near their centers and 
approach p{r) oc only beyond a few thousand AU 
(WSHA, AWM96, Fig. 1). In some cases, sharp outer 
boundaries, much steeper than p{r) oc r“^, are ob¬ 
served at a finite core radius (e.g. Abergel et al. 1996). 
This raises the possibility that, in some instances at 
least, the initial conditions for protostellar collapse de¬ 
part significantly from a singular isothermal sphere. 

2. The outflow momentum flux of embedded protostellar 
sources correlates very well with their circumstellar en¬ 
velope mass (see Fig. 2 and Fig. 7). This suggests that 
the mass ejection and mass accretion rates of proto¬ 
stars both decline with time during protostellar evolu¬ 
tion from Class 0 to Class I sources (BATC). 

3. Recent self-consistent hydrodynamical calculations of 
protostellar collapse indicate that initial conditions 
characterized by centrally flattened density profiles 
lead to a transitory phase of enhanced accretion im¬ 
mediately following the formation of the central hy¬ 
drostatic protostar. This behavior is found regardless 
of whether the influence of magnetic fields is ignored 
(FC93) or fully accounted for (Tomisaka 1996). The 
nature of the transitory accretion peak is, however, 
left somewhat unclear by these numerical simulations. 

4. On the basis of (3), we propose that (2) is a direct 
consequence of (1). 

5. In order to elucidate the physical origin of (3), we 
analytically follow the history of protostellar accre¬ 
tion, using Lagrangian calculations based on a new 
type of self-similar ‘gravity-dominated’ collapse solu¬ 
tions (Sect. 3.2, Sect. 3.3, and Appendix A). We claim 


that these new solutions provide a better description of 
the collapse than the usual ‘pressure-dominated’ self¬ 
similar solutions (e.g. Shu77, WS85), when pressure 
becomes negligible in the self-gravitating flow some¬ 
time after the onset of collapse, prior to stellar core 
formation. Comparison with the numerical simulations 
shows that the gravity-dominated solutions are ade¬ 
quate, at least qualitatively, in the supersonic region. 

6 . In agreement with (1), we start our supersonic cal¬ 

culations at < = to < 0 from an idealized pre-stellar 
core consisting of a strictly flat inner plateau up to a 
radius r^r, an ‘envelope’ up to r^, and a steeper 
power-law ‘environment’ farther out (Fig. 3). In our 
gravity-dominated description, the central plateau re¬ 
gion first collapses homologously to form a finite-mass 
hydrostatic stellar core at t = 0. Observationally, this 
initial phase, which does not exist in the standard Shu 
picture, should correspond to ‘isothermal protostars’, 
i.e., supersonically collapsing cloud fragments with no 
central YSOs. It is followed at t > 0 by the main accre¬ 
tion/ejection phase, during which the non-zero central 
point mass accretes the surrounding envelope. As long 
as the gravitational influence of the initial point mass 
is significant, the accretion rate remains higher than 
the Shu value, It then quickly converges to¬ 

wards aljjfG. At late times, accretion of the outer 
environment leads to a terminal phase of residual ac¬ 
cretion/ejection, during which the accretion rate de¬ 
clines below the Shu value (see Fig. 5). 

7. We use our analytical model to fit the diagram of out¬ 
flow efficiency Fco versus normalized envelope 

mass obtained by BATC (see Fig. 7). To 

the extent that there is a direct proportionality be¬ 
tween accretion and ejection, this diagram should pro¬ 
vide an empirical measure of the accretion history of 
the sampled protostellar objects. A good overall fit is 
found when the model boundary radius rf,, is not much 
larger than the radius tn of the flat inner plateau, i.e., 
rf, ~ 1.6rjv- This requires that the fraction of cloud 
mass in the central plateau region be relatively large, 
MN/Mcioud ~ 30 %. (However, since the collapse is in 
fact less violent than in our pressure-less calculations, 
these values of rf, and Mn should be taken as indica¬ 
tive, and are likely to represent only a lower limit and 
an upper limit, respectively.) 

8 . Based on (7), we tentatively associate the short pe¬ 
riod of energetic accretion/ejection predicted by our 
model at the beginning of the accretion phase with 
the observationally-defined Class 0 stage (Sect. 4.1 and 
Fig. 7). In this view. Class I objects are more evolved 
and correspond to the longer period of moderate accre¬ 
tion/ejection when the accretion rate approaches the 
Shu value. (See, however, point 10 below.) 

9. We also find that the Ophiuchus F’coc/Tboi versus 
Menvl^ diagram differs markedly from the cor¬ 
responding diagram in Taurus: a clear contrast be- 
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tween Class 0 and Class I objects is observed in Ophi- 
uchus (Fig. 9) which is not seen in Taurus (Fig. 10). 
This points to different accretion histories in these two 
nearby star-forming clouds, which we interpret as aris¬ 
ing from differences in initial conditions (Sect. 4.3). 
Both outflow and dense core observations suggest that 
the relative importance of the central plateau region 
in the initial density profile is significantly larger in 
Ophiuchus {tn ~ rb, Mm/M doud ^ 30 % according to 
our approximate model) than in Taurus (r^r ~ 0.2rj,, 
Mm/M cioud < 10 %). 

10. According to our model fit of the Ophiuchus Fco c/Lboi 
versus Me„v/Tboi°'® diagram, most of the p Oph Class I 
YSOs should be in their terminal accretion phase (see 
Fig. 9). This is not surprising since the relatively small 
fragmentation length scale observed in p Oph implies 
that only a finite reservoir of mass is effectively avail¬ 
able for the formation of any given protostar. 

11. In conclusion, the ‘standard’ theory of Shu and co¬ 
workers appears to describe protostellar evolution 
quite satisfactorily in regions of isolated star formation 
like Taurus. In our view, this is because in these regions 
most stars probably form following the self-initiated 
contraction/collapse of dense cores due to ambipolar 
diffusion. 

However, in star-forming clusters such as p Ophiuchi, 
the standard Shu theory is less appopriate since it does 
not account for fragmentation and multiple star for¬ 
mation. In these regions, star formation may be in¬ 
duced by the impact of (slow) shock waves (e.g. Boss 
1995), and protostellar cores may form by supersonic 
(or superalfvenic) implosion of dense clumps, rather 
than slow (subsonic) ambipolar evolution. In this case, 
the collapse/accretion history advocated in the present 
paper (Sect. 3.3) and summarized in (6) is likely to 
provide a better description (see Fig. 9 and Sect. 4.3). 
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Appendix A: Lagrangian collapse and self¬ 
similarity 

In Henriksen (1989 - H89) it was shown in effect that in 
terms of a Lagrangian fabef r defined as the radial position of 
an isothermal spherical shell at a hducial instant t — to when 
the sphere is at rest, and a variable 

^ = y/2GM(r)/r3 {t - U) (Al) 

{M{r) is the mass initially inside the radius r), the subsequent 
position of the shell R{t,r) can be written as 

R = rS{i,r). 


The ‘scale factor’ S obeys the equation 

(de5)^-l/5 + a.(r)^H(C,r) = -l. (A2) 

Here, we have dehned = afr/GM{r), where as is the 
sound speed, and Ti = \a p{^,r)/po{r). The function po{r) is 
the initial density (at t = to) that is related to M(r) through 
Po = ( 47rV^ ) shall discuss the possible forms of this 

function further below. 

In order to obtain a complete description of collapse from 
rest, we must add the equation for the density in the form 




SM1 + 
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dlnS __ 

51nr -D(r) d 


InjSA 

Ine J 


(A3) 


where the function .^ = ||^ = § - l). 

These equations are equivalent to the equations used pre¬ 
viously in Eulerian form by most other authors (e.g. FC93). 
Clearly S{f = 0,r) = 1 is one boundary condition and 
(5 = 0, r) = 0 is another, but in addition we must impose an 
equilibrium condition if this is to be a ‘natural’ starting point. 
The equilibrium condition consists in requiring Po{r) to satisfy 
the Lane-Emden equation for an isothermal gas sphere. Con¬ 
sequently it is either the general function po/po{0) = 
discussed for example in Chandrasekhar (1939) and used in 
the Bonnor-Ebert discussion of instability (e.g. Bonnor 1956), 
or it is Po = X/r^ everywhere, which is the singular isothermal 
sphere (SIS). In the latter event = 1/2 is constant and 


A — cTd 


AtvG’ 


(A4) 


Moreover the variable ^ is recognizable in this case as 


(H5) 

The SIS is exceptional in that Eqs. (A2) and (A3) show 
that the development can be immediately ‘self-similar’; that 
is <S(^,r) = S{^), and Ti. = H{^) since D{r) = 1 for all r. 
This is the initial condition and collapse hrst studied by Shu77 
and subsequently developed by Shu and co-workers into the 
‘standard model’ of protostar formation. 

However, a closer inspection of Eqs. (A2) and (A3) com¬ 
bined with our knowledge of the expected parametric collapse 
instability discovered by Bonnor and Ebert reveals that there 
are several other possible lines of development. These are all 
based on the assumption that the initial equilibrium is not 
that of the SIS, but rather a classical Bonnor-Ebert solution 
with a flattened core (po(0) is finite and ip ~ x^/6 for small 
X = 47rGpo(0)/air) and a r~^ halo. 

In general, if we ignore the somewhat artihcial means of 
launching the collapse by a sudden cooling of the gas so that 
Og ^ 0, the development from the Bonnor-Ebert sphere is 
clearly not self-similar since as and D are complicated func¬ 
tions or r. But our fundamental theoretical remark in this pa¬ 
per is that a self-similarity in terms of the general variable 
i{t,r) (as in (Al)) can develop dynamically if ^ 0 as the 
collapse develops. For then S S{^) by Eq. (A2) and the den¬ 
sity tends to p/po = S'”®/)! — -p^dlnS/dln^) as used in the 
text. 

To place this idea in context it is worth remarking that none 
of the entire ‘hyperbolic’ (that is characterized by a wave prop¬ 
agation speed) family of self-similar solutions found by WS85 
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can begin exactly from the Bonnor-Ebert initial condition ei¬ 
ther. They must then arise dynamically in the same way that 
we suggest occurs for the zero pressure flow. The ‘band zero’ 
solutions may represent an intermediate stage however because 
of their flattened, dis-equilibrium, central peaks, as discussed 
in the text. Such solutions are all constrained to pass through 
an instant where p oc r~^ because of their assumption of the 
strict self-similar symmetry (see e.g. Carter & Henriksen 1991) 
expressed in Eq. (A5) (normally authors set to — 0). This in¬ 
stant occurs at ^ = 0 (see e.g. Eq. 1 of Hunter 1977) which is 
normally taken to coincide with t = 0. This family of solutions 
may be discussed elegantly in terms of the Lagrangian formu¬ 
lation as discussed in H89. However since these solutions are 
well known (WS85), we pursue the novel idea in this paper as 
noted above, that the self-similar symmetry itself may change 
dynamically from that of (A5) to that determined by ^ where a 
more general mass or density profile is allowed. Essentially we 
suppose that the symmetry changes from being ‘hyperbolic’ 
to ‘elliptic’, which is the symmetry appropriate to the local 
acceleration that is present in a spherical flow dominated by 
gravity. 

This requires that the flow evolve away from the initial 
equilibrium under the influence of a Bonnor-Ebert-Tomisaka 
type instability towards a density profile where is small 
compared to — 1 (Eq. A2), at least in some hnite region. At 
this instant, which we designate as —tff{N) in the text, where 
the gas will no longer be in equilibrium but rather in supersonic 
flow {{dS/d^)'^ » Os), the current position of the mass shells 
becomes the new comoving variable r. Moreover S S{^) 
becomes determined by the pressure-free version of Eq. (1) in 
the text, for which we take the bound solution that becomes 
Eq. (4) of the text near ^ = tt/2. Here we distinguish between 
to and —tff (N) to allow for an uncertain duration between the 
actual initiation of collapse and the onset of supersonic flow (in 
the text the formulae are simplified by taking to = —tff{N)). 

This scenario will be realized in the vicinity of rjv for ex¬ 
ample if at —tff{N) the density profile is roughly as in Fig. 3 
of the text, and if a^VN/GM n << 1- We assume that this 
condition is attained during the evolution between t = to and 
t = —tff{N). From this point on we describe the flow as pos¬ 
sessing the elliptic self-similar symmetry of zero-pressure self- 
gravitating flow. 

One may note that the L/P solution, which seems to be 
a better (although not perfect) fit to the simulation results 
than is the S solution (FC93), also develops with as small. In 
fact such a band 0 solution (WS85) is described by Eqs. (A2) 
(with k = 0), (A3) and (A4), and so we may calculate the 
radial velocity near the singularity as 


Appendix B: Asymptotic value of the accretion rate 


We may make contact with the accretion rates found by 
Shu77 and by FC93 by modifying slightly the argument of H94. 
If as above rjv is the initial radius at which the profile becomes 
sufficiently flat to initiate the homologous nuclear collapse then 

7, r 47r 3 

Mn = —PNrN- 

Moreover we dehne here the equivalent of the ‘virial’ radius 
used in H94 by referring to the dimensionless variable (j" = 
r^/AirGpcIas. This gives in the present context 


rjv = Cjv 


tts 

7,/A'kGpn 


(Bl) 


If the initial sphere were exactly Bonnor-Ebert, then the differ¬ 
ence of density from pjv at ^jv is oc (Cv/6 (e-g- Chandrasekhar 
1939). Thus Cjv < 1 to satisfy the flatness criterion. However 
we have detected a tendency for the density profile to flatten 
on the outside so it is possible that ^ > 1 at the initiation 
of homologous collapse. Thus we retain C,m here as a param¬ 
eter. As remarked in H94, the preceding two relations can be 
rearranged into the useful forms 


, 3 , GMn 
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and 
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ATTr% 
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The units in terms of which we have calculated our accretion 
rate now follow as 


Mn = 4nrNpN \J 2 GMn / rjv 


2^3 at 
G’ 


(B4) 


Since moreover the asymptotic value of the accretion rate in 
these units was found above to be — \/3 (2/7r), we predict for 
the accretion plateau that is established as the r~^ outer region 
falls in the rate 

= (55) 

TT (_T 

The numerical factor is 0.900 Cv so that ((jv = 1.027 in order 
to yield the asymptotic accretion rate predicted by the Shu 
solution (a numerical factor of 0.975 in Eq. B5). This value of 
C,N seems quite possible in view of the expected flattening. 
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